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ABSTRACT 


Fueled by ongoing rapid decline of biodiversity (T), ecology is in the midst of a lively debate on the effect of 
species loss or introduction on food web stability |2l[3). In food webs, many interacting species coexist despite 
the restrictions imposed by the competitive exclusion principle 0H) an< I apparent competition (6). For the 
generalized Lotka-Volterra equations Q, sustainable coexistence necessitates nonzero determinant of the 
interaction matrix. Flere we show that this requirement is equivalent to demanding that each species be part 
of a non-overlapping pairing, which substantially constrains the food web structure. We demonstrate that a 
stable food web can always be obtained if a non-overlapping pairing exists. If it does not, the matrix rank can 
be used to quantify the lack of niches, corresponding to unpaired species. For the species richness at each 
trophic level, we derive the food web assembly rules, which specify sustainable combinations. In neighboring 
levels, these rules allow the higher level to avert competitive exclusion at the lower, thereby incorporating 
apparent competition. In agreement with data HH9J, the assembly rules predict high species numbers at 
intermediate levels and thinning at the top and bottom. Using comprehensive food web data QOllIIL we 
demonstrate how omnivores or parasites with hosts at multiple trophic levels can loosen the constraints and 
help obtain coexistence in food webs. Hence, omnivory may be the glue that keeps communities intact even 
under extinction or ecological release of species. 
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INTRODUCTION 


In food webs, complexity arises from combining a large number of species (the nodes) and a large number of 
relations between these species (the links). Addressing the latter, recent attention was devoted to the structure 
of links using e.g. the random, cascade and niche models lT2l stirring a prolific debate on the role of the 
link distribution regarding food web stability meannua. We take a complementary approach: Using standard 
consumer-resource equations, we demonstrate fundamental constraints on node diversity in a food web, termed 
food web assembly rules. 

For consumer-resource relationships, the competitive exclusion principle states that when two consumers com¬ 
pete for the exact same resource within an environment, one consumer will eventually outcompete and displace 
the other (4j (5| . It is known that the number of coexisting species cannot exceed the number of resources these 
species compete for m. Expressed more generally, the number of coexisting species cannot be greater than the 
number of distinct regulating factors in the community E3- For trophic communities of several levels, it was 
subsequently stated that the number of species on any trophic level could not exceed the sum of the numbers on 
adjacent levels mi- Experimental studies do demonstrate strong correlations between consumer and resource 
diversity UU [ 20 ] EUI 22 ) . These observations highlight that also the consumer plays a critical role in shaping the 
network of species, even when direct interaction between resource species is absent, an observation captured in 
Holt’s paradigm of apparent competition |6). 

Despite the existing theoretical constraints and empirical findings, a selective theory for stable coexistence 
of many species in food webs is currently lacking. This lack may partially be due to the complexity of the 
many-species interactions, yielding an uncontrollable number of parameters and hampering direct calculations 
or simulations of sufficient generality. Notwithstanding these complications, progress can be made when nec¬ 
essary conditions are demanded. For an ecology, consisting only of a resource and a consumer level, we have 
recently shown that coexistence requires that the species richness of both levels is balanced and that a cascade of 
parameter values must be maintained ll23l . Examples of such systems may be the phage-bacteria ecology in the 
Atlantic Ocean ll24l[25 i or laboratory ecologies. 
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However, in food webs, a subset of trophic levels can generally not be considered in isolation. A species’ niche 
is determined by its entire set of interactions, which generally may be composed of both beneficial and harmful 
interactions, i.e. the species may act both as a consumer or resource. Further, many food webs contain omnivorous 
interactions, i.e. those where one species preys on several other species that are located at more than one trophic 
level. To derive necessary conditions for coexistence in food webs, a more general starting point is required. 

Based on the generalized Lotka-Volterra equation 0, we here show that in sustainable food webs each species 
must be part of a non-overlapping pairing. We define a non-overlapping pairing as a topological pattern for a 
directed network, where each species must contribute to a closed loop and none of the loops may overlap (Details: 
Methods). Mathematically, this is a consequence of demanding that the determinant of the interaction matrix be 
nonzero. We then analyze the implications of this pairing for the species richness at different trophic levels. We 
make predictions for secondary extinctions and assess the stabilizing effect of parasitism and omnivory. 

METHODS 

Steady state equations. For consumer-resource interactions in food webs, the generalized Lotka-Volterra equa¬ 
tions 0 are 
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for species at trophic levels l > 1. We distinguish a species by the set of links that connect it to predators and 
prey or nutrients and the strength of these links (Details: Sec. S4). In Eqs[I]and|2] S, <li with i = 1,..., ni are the 
densities of species residing at trophic level l, m is the species richness at level l, denote the maximal growth 
rates of pji describe differential consumption of the basic resources by the S^\ denote the decay rate of 
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species S^\ r/^f 11 are the interaction coefficients between a species S' 1 / 1 with species Sf 1 \ and ft 1 / 1 are the 
efficiencies of reproduction of species S^ 1 when consuming species Sm 


In the steady state, the time derivatives S^ and S'® on the LHS of Eq.Q] respectively Eq.[2] vanish and we 
have the equations 



ki~ \ respectively 


(3) 
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Collecting all constant coefficients (RHS of Eqs[3]and[4} in the vector k and all interaction coefficients on the 


LHS in the interaction matrix 1Z, we have the linear matrix equation TZ ■ S = k, where S is the vector of all species 
densities. For completely shared nutrients, the competition factors pj, = 1. (Details: App. Secs SI and S2). 
Parasite interactions. Parasites have complex life-cycles that can demand several hosts (2611271. Notwithstanding 
these complications, we here formally treat them as consumers, respectively resources of independently acting 
other species. Also, we simplify concomitant links in terms of simple linear responses (details: SI Sec. S9.3.3). 
Non-overlapping pairing. det(JZ) f 0 can be fulfilled if every species is paired with another species or nutrient 
(this constitutes a perfect matching (28)). For food webs with sharp trophic levels, the network is bipartite and 
therefore it is required that such a perfect matching exists. When species with variation in food chain length are 
present, one may generally obtain nonzero det(JZ) by covering the entire network with dosed loops of directed 
pairings (i.e. cycles). This is a sequence of nonzero matrix elements IZij, TZjk, 7 Zki, ■ ■ ■, 1Z m , (Sec. S6), i.e. a 
chain of directed pairings where the direction is maintained and the last element connects to the first. A directed 
pairing represents one nonzero matrix element, whereas a pairing also includes the symmetric element. 

Empirical food web data. We use high-resolution data on seven food webs including free-living and parasite 
species: The North American Pacific Coast webs Carpinteria Salt Marsh (CA), Estero de Punta Banda (PB), Bahia 
Falsa (BF) |29]|T0j; the coastal webs Flensburg Fjord (FF) lf30l . Sylt Tidal Basin (ST) 071 . and Otago Harbor 
(OH), New Zealand l32l . as well as the Ythan Estuary (YT), Scotland 11261 . These food webs describe consumer- 
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resource interactions between basal, predatory and parasite species. A compilation of all seven food webs has 
recently been provided mm. Specifically, the data distinguish three types of links: (i) links between free-living 
species only (“Free”), (ii) additional links between parasites and other species (“Par”) and (iii) links from free- 
living consumers to the parasites of their resources (“ParCon”), i.e. so-called concomitant links. Details on data 
analysis: SI Section S9. For the empirical data, the lack of niches, i.e. nullity d = S — rank(7?,), was computed 
by using random values for all nonzero entries of the respective matrix 1Z. Basal species were each given an 
individual nutrient source. In “ParCon sym” a subset (20 percent) of concomitant links were randomly selected 
to be symmetric (Details: Sec. S9.3.3). In the data analysis and simulations (Fig. |5] and |6j, the trophic level of a 
species is defined by its prey-averaged food chain length (Sec. S9.2.1). 

Simulations. We perform two types of simulations: (i) In-silico assembly of a tree-like food web (Fig. [3}:), where 
parameters are chosen according to constraints discussed in SI Sec. S2. The numerical values of the parameters 
are: For the interaction and growth coefficients r) = /? = k = 1 for all links present (solid black arrows in Fig.[T|y), 
as well as the decay coefficients {ai,..., c^s} = {.1, .1, .16, .1, .12, .15, .1, .1}, where the labels are as indicated 
in Fig. [ 3 J 3 . Each new species is introduced at low density and time-integration is continued until steady state is 
reached (using Mathematica NDSolve method), (ii) An idealized food web was constructed by using the average 
species counts at levels n ? ; obtained from all empirical data sets and initially assuming sharp trophic levels for 
all species. Sharp trophic levels were obtained by rounding each species’ chain length to the nearest integer value 
(Fig.[6ji). With the constraint of these trophic levels, a number of links was assigned to match the empirical average 
for free-living food webs (Fig. [6}t). When adding further species, the empirical average of parasite species count 
was used (47 species). To obtain Fig. [6] initially, each parasite received one link. In the cases (c) and (d), this link 
connected the parasite to any existing free-living species. In the cases (e) and (f), this link connected the parasite 
to any existing species at trophic levels 3 or 4. For any subsequent link, a parasite was chosen at random. A link 
was then formed in three ways: Case c: randomly to connect with another existing species; Case d: randomly, but 
only to species at the same trophic level as for the initial link of that parasite; Case e: randomly to other existing 
species at levels 3, level 4; Case f: randomly to other existing species at levels 3, level 4 or another parasite (More 
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detail: Sec. S10). 


RESULTS 

Theory 

We describe the interaction of species on L trophic levels by the generalized Lotka-Volterra equations IT71 l34ll . 
Basal species are constrained by the system carrying capacity while the consumers are assumed not self-limiting, 
and trophic interactions occur through the linear type-I functional response (Methods). Such equations have been 
widely used in community assembly models, where food web networks are assembled by numerically analyzing 
the equations to find parameter sets with stable and/or permanent coexistence solutions |35l [36 1 371. Here we 
take an alternative path by first finding a necessary condition for coexistence in terms of species richness, which 
results in the food web assembly rules that constrains the network topology. We subsequently show that when 
these assembly rules are fulfilled, stable network structures can always be obtained. 

In the steady state, we have the matrix equation 1Z • S = k, where S is the vector of all species densities, 7Z 
is the interaction matrix between the species, and k is the vector of growth and decay coefficients. Note that 7 Z 
has a block structure with nonzero entries only for interactions between neighboring trophic levels but not within 
the same level, and that the positions of these matrix elements are symmetric due to mutual interaction between 
predator and prey (Fig.jT^). Stable/permanent coexistence requires that a feasible solution S* = 7Z 1 k > 0 exists 
03J [36]]. Structural stability, i.e. robustness against parameter perturbations, of a feasible coexistence solution 
requires that det(TZ) be nonzero (Sec. S2), a condition required for the existence of the matrix inverse in a linear 
equation l38l . In the following, we analyze what this basic condition means for the species richnesses at the 
different trophic levels of a food web. We then discuss also the stability of a given solution. 

We first specialize to the case of a single shared basic nutrient, e.g. sunlight, used by all primary producers 
s [ 11 , with 7ii the species richness on the first trophic level and i ranging from 1 to n-j . This corresponds to setting 
the competition terms pji = 1 for all 1 < i,j < n\ in Eq. [Q thereby yielding a block of ri\ x v,\ unit entries 
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in the lower right block. (Example: Sec. S4). Nonzero determinant is achieved if it is possible to identify a path 
of matrix elements that only contains elements from the non-zero sub-matrices bordering the diagonal (Fig. [Qo. 
This is equivalent to demanding that every species be part of a consumer-resource pair connecting neighboring 
trophic levels and none of these pairs overlap (known as perfect matching of a bipartite graph in graph theory 
l28l ). The pairing guarantees that no species share exactly the same niche, i.e. a particular set of interactions 
with resources and consumers (Sec. S4), and manifests the competitive exclusion principle (Fig. [JJ?, i and ii). For 
primary producers, pairings may involve the nutrient source (Fig. nk inset). In that case, at least ni — 1 species 
at level two are required for pairing of the remaining basal species. We term this structure resource-limited. In 
this case, ri 2 > «i — 1 is required. In the example (Fig. [TJt) this condition can indeed be fulfilled, because a 
sufficient number of species exists at level two. If pairing with nutrients is not used, ri 2 > n \ is required — a 
consumer-limited configuration due to the biomass restriction imposed by consumer predation (Fig. S2). In turn, 
those species left unpaired at level two must be paired by species at level three. 

Defining N 0 and N e as the sums of node richness at odd and even levels, respectively, it is then easy to check 
that the general constraints become 

A = N 0 -N e G {0,1} , (5) 

which encompasses the competitive exclusion principle IU[5](23)- Defining L as the top trophic level, for each of 
the two options in Eq.[5]a set of L — 1 nested inequalities arises relating the species counts n % : 


(<) 

Til > 
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{ii) 

ri2 > 

ni — A 
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n 3 > 

ri2 — n\ + A 

(6) 

O) 
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IV 

n 3 — n 2 + ni — A 



This sequence continues until nz, is reached. The case A = 0 is consumer-limited, whereas A = 1 signals resource 
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limitation. Note that these rules allow simple assessment on which food webs can have coexistence-solutions. To 
illustrate this, we give two simple examples of food webs that cannot coexist by Eas l5l6l (Fig.lTb. iii and iv). Our 
rules thereby are more selective than those in previous work fl8l . There, a requirement was stated for trophic 
communities of several levels, where the number of species on any level could not exceed the sum of the numbers 
on adjacent levels — a condition that would e.g. not rule out the two simple webs (Fig. []J>, iii and iv), hence not 
detect the lack of niches. 

The conditions (Eqs l5l6b rule out stable coexistence when violated, yet they leave unanswered if all populations 
can be positive or stable when they are met. For any choice of species counts fulfilling Eas 15161 a feasible solution 
can be obtained by starting from the complete non-overlapping pairing of species, i.e. a subset of links. A tree-like 
backbone of links will be obtained, which preserves the complete non-overlapping pairing. Additional links give 
each species access to nutrients. For this structure, we find that the parameters can always be assigned to yield 
positive steady-state population densities for all species. This is possible by enabling trade-offs between the invader 
and the resistant species, so that the fitness of any invader, e.g. described by its decay coefficient, is sufficiently 
large to maintain competitiveness of the residents lf39i . An example is shown in the subsequent section. 

The existence of a Lyapunov function is a sufficient criterion for global stability [40]. Using this, in extension 
of the analysis of simple chains of species ED, we demonstrate that these food webs are always stable, i.e. 
the Lyapunov function decreases monotonically with time after the initial perturbation (Sec. S2.2). We thereby 
establish that for each combination of species richnesses consistent with the conditions (Eqs 1516b . a feasible and 
stable food web exists. 

In practice, there may be several basic nutrients, such as different chemical compounds or spatial or temporal 
subdivision fl42l . If ns > 1 separate nutrients are available, the assembly rules yield ns +1 sets of conditions 
analogous to those in Eq. [6] where any A E {(), ns} is allowed (Sec. S7). The presented theory assumes simple 
predator-prey couplings. Non-linear interactions (i.e. type-II response) and cannibalism can be included by adding 
diagonal matrix elements in Fig. [TJi, corresponding to species that pair with themselves (Secs S7,S8). Another 
future extension of the theory could be that of incorporating frequency-dependent predation |43l . 
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Food web assembly 


We return to a single nutrient (ns = 1) and now discuss species richness at the different trophic levels. In the 
consumer-limited case, Eq. EJi) restricts n\ to an upper bound given by 112 . Eq. [6jzz) conversely restricts n -2 to 
numbers equal to the sum of the total species richness in both neighboring trophic levels, hence the limitation to 
the abundance of n ,2 is much weaker compared to that of ri \. Accordingly, the basal level cannot constitute the 
global maximum of node richness within a level. Moving further, ns could again exceed 11.2 but possible limitations 
stemming from the count of top predator species become noticeable. When starting at the top predator level L, 
by symmetry, similar constraints as in Eq. [6] hold: ul — n-L-i < 0. Together with Eq. [6jz) we have for the 
consumer-limited state, ni+nL < 112 + n l—i, i.e. species at intermediate trophic levels generally dominate food 
web biodiversity (Fig. 0. The condition is similar for the resource-limited state, where the limit is shifted by one 
species. For increasing L, these equations predict further increase in the contribution of intermediate species to 
total biodiversity. 

Given these general constraints, we can now discuss food web assembly. Consider a simple food web and its 
interaction matrix (Fig. [3b). Graphically, an allowed structure is again manifested by permitting non-overlapping 
pairing of species. Food web growth is characterized by alternating transitions between resource- and consumer- 
limited states. Initial community growth requires the presence of a single primary producer (Fig.[T|“>). hence n 1 = 1. 
The only possible addition is then a consumer, preying on the existing producer. This entails immediate addition of 
a trophic level (m = ri 2 = 1) and transition to a consumer-limited state. The assembly rules subsequently permit 
additions at different trophic levels, but alternation between consumer and resource limitation must be preserved. 

For each food web consistent with the conditions in Eqs 15161 parameters can be assigned in a way that food 
web evolution is possible, i.e. that feasibility and stability are even achieved after each addition of species. To 
exemplify this, we assign parameters for all eight species of the example, tree-like, food web (Fig. [ 3 J 5 , |[39l ). 
Using in-silico simulations, we mimic the evolution of this web, starting from only a single species (Fig. [3};). Any 
addition of new species leads to a transient disruption of population densities, but a new feasible steady state is 
eventually established (Fig. [3}:). To quantify stability, we further compute a Lyapunov function, and verify that it 


10 


monotonically decreases after each addition within each stage of the food web (Sec. S2). We have checked that 
addition of further weak links is compatible with stability and feasibility, hence in principle allowing arbitrary link 
structure. 

Our assembly rules are easily generalized to food webs containing omnivores or parasites with hosts at several 
trophic levels (Sec. S6). Consider again Fig. [3jt, but imagine that another species is added on trophic level 4, 
causing a violation of the assembly rules. This violation can be rectified by a single generalist omnivore (Fig. [4]), 
corresponding to an additional row and column of nonzero coefficients in the interaction matrix. The omnivore can 
be interpreted as a consumer preying on all trophic levels. Hence, one can choose a level i and let its species count 
rii be increased by one unit to again satisfy the assembly rules. In general, omnivory or parasitism can make the 
graph non-bipartite. If so, the non-overlapping pairing to achieve det(lZ) ^ 0 is extended to covering the entire 
network with closed loops of directed pairings (Sec. S6). 

Data analysis and simulations 

What do the assembly rules teach us about real food webs? For seven detailed empirical food webs (Details: 
Materials and Methods) containing both free-living and parasite species, we determine the difference between the 
respective total number of species S and the maximum number of linearly independent rows, d = S — rank(T^) 
(Details: Methods). Linear dependence can be seen as the sharing of a specific niche by several species, therefore 
d measures the lack of niches in the given food web. When d = 0, the assembly rules are satisfied ( det(JZ ) ^ 0). 
As done previously cm we distinguish webs formed by: links between free-living species only; links between all 
free-living and parasite species; with additional concomitant links (Methods and SI Sec. S9.3.3). 

The free-living webs have substantial structure, an example is given in Fig. [5ji and others in Sec. S9. Many 
species are located at sharp trophic levels (Fig. S7), a feature that manifests itself by the blocks of white spaces, 
i.e. absent interactions (Fig. [5}i). The organization into sharp trophic levels entails stricter demands on the com¬ 
binations of species richnesses (Eas 1516b . and many sub-webs consisting only of free-living species do not seem 
to fulfill these demands. This is quantified by an associated lack of niches in most free-living webs, i.e. d > 0. 
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Additional analysis reveals that all free-living webs are in a consumer-limited state, i.e. species richness in even 
trophic levels dominates (Sec. S9). 

We contrast these findings with simple models, namely the cascade fl2l and niche models lfT3l . Using number 
of species and links from the empirical datasets, we generate network samples (Fig.[5j),c). The resulting interaction 
matrices are characterized by very little structure in terms of trophic levels (the white blocks are all but missing). 
When repeating the simulations for all seven webs (Sec. S9) and obtaining the corresponding rank deficiencies, we 
find that the networks simulated using the cascade and niche models consistently give dfs 0 and are less structured 
than empirical data (Sec. S9). 

We further quantify the organization of species by the chain length distributions for the empirical and modeled 
networks, where much broader distributions are found for the models. Quantifying omnivory by the standard 
deviation of chain lengths ll44l . modeled networks consistently yield substantially higher fractions of omnivory. 

We now consider the webs involving parasites (Fig. |5}i,e). At the edge of the panels we indicate by a color¬ 
coding, where, in terms of trophic level, parasites enter and how the free-living species are re-organized. Notably, 
parasites predominantly enter at high trophic levels, the lower section (approximately levels one and two, red to 
green colors) remains nearly unaltered by the inclusion of parasites. Specifically, interactions of a given parasite 
generally involve several free-living species at multiple trophic levels (Sec. S9), acting to loosen the structure at the 
higher levels ( 17,3 and 714 , compare Fig.[5ji,d) and reaching agreement with the assembly rules (Sec. S 6 ). In other 
words, some of the parasites can be seen as effectively acting as odd-level species, thereby relaxing the initially 
consumer-dominated free-living webs to a more balanced state. 

Concomitant links (Fig.|5^) cause further entanglement of trophic levels, open additional options for possible 
pairings and systematically increase sustainability of food webs. Concomitant links require additional consider¬ 
ation, as they generally are directed links where a parasite is consumed by its host’s predator, i.e. they denote 
a detrimental effect on the parasite population. A positive impact on the predator’s population may however not 
always result. We call such links asymmetric concomitant links. If the predator’s population does benefit from the 
consumption, such as observed in some studies Il44l f45l. we use the term symmetric concomitant link. Nonetheless, 
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such directed links can lead to additional non-overlapping pairings, when a closed loop of directed links is formed, 
e.g. a triangle (loop of length three). In the empirical webs it is noticeable that inclusion of asymmetric con¬ 
comitant links only rarely yields rank improvement. Investigating this further, we find that many food webs with 
parasite interactions already contain sufficient numbers of loops to allow pairings between parasites and free-living 
species. The limitation arises because a surplus of parasites exists. Each loop will involve at least two free-living 
species but only one parasite, making it impossible to find non-overlapping pairings for all parasites. For those 
webs, only the inclusion of symmetric concomitant links leads to an additional improvement of rank, since then 
each parasite can be paired with a single free-living species. In two webs, where loops are rare, even asymmetric 
concomitant links help improve the rank. Details and simulations are available in SI sections S9.3.3 and S10.4. 

For all empirical food webs, we summarize the effect of the different link additions on rank deficiency d 
(Fig. |3). We find a general decrease of d as more parasite links are added. Notably, full rank is sometimes not 
achieved, even when all available concomitant links are included. Indeed, empirical food web datasets may often 
be incomplete, as some links can be difficult to detect. Our findings may serve as a means of identifying possibly 
missing data, most notably in the Ythan Estuary food web, where overall link density is low and parasite-parasite 
links are completely absent (further details on individual food webs: Sec. S9.3.2 and S10.3). 

As mentioned above, existing food web models generally produce network structures lacking rank deficiency 
(d « 0), even for the free-living webs. When instead starting from model networks with similar link density as 
the empirical webs but sharp trophic levels, we obtain d substantially larger than zero (Fig. 13- Parasites, with 
their complex life-cycles ll26l . often consume species at varying trophic levels during different stages of their lives. 
Adding species that each interact with species at multiple trophic levels (Fig. [3), he. mimicking the addition 
of parasites, our modeled webs give systematic decrease of d as links are added (black line in Fig. (3)- When 
adding species that each interact with a single (Fig. [3). or exclusively higher trophic levels (Fig. [6^), saturation at 
d > 0 occurs (blue and red lines in Fig.|3>)- When simulating also parasite-parasite interaction, d is also found to 
decrease (Fig. [3 Simulation details: Methods and Sec. S10). Overall, these simulations suggest that addition of 
random species without the feature of interaction with multiple trophic levels is not sufficient to explain removal 
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of rank deficiency. 

We have further explored addition of omnivorous links to the free-living web. The rank deficiency d is reduced 
rapidly if the addition happens randomly for all trophic levels, but the reduction is limited if omnivorous links 
occur mostly at the trophic level 3 as in the real data (Sec. S10.2). We have also performed extensive simulations 
on the effect of concomitant predation, which further emphasizes the importance of parasite-parasite interaction in 
achieving coexistence for some webs (Sec. S10.4). 

DISCUSSION 

The food web assembly rules generalize the competitive exclusion principle to food webs of any number of species 
and trophic levels. They quantify which combinations of species richness at the different trophic levels can yield 
coexistence solutions. We show that for any of these combinations, there are stable and feasible network structures. 
Demanding full rank of the food web interaction matrix expresses the simple notion that each species must occupy 
a unique niche and leads to biologically plausible combinations of species richnesses at the different trophic levels. 
The requirement, i.e. non-overlapping pairing or equivalently nonzero determinant, is simple and directly allows 
us to evaluate the self-consistency of empirical data. The rules help explain that actual food web networks are far 
from random and more structured than those obtained from the traditional niche and cascade models. While some 
food web datasets do fulfill our conditions, networks known to lack interactions, e.g. the Ythan Estuary web j26l . 
stand out as particularly far from reaching agreement with the rules. 

One immediate consequence for food webs with species predominantly organized according to trophic levels, 
e.g. many free-living webs, is that species richness at the basal and top-predator levels should be limited by the 
species richness of the respective neighboring levels (compare Eq. [6j. This can explain the observations, both for 
terrestrial 0 and marine food webs 0 , which report greatest species richness at intermediate trophic levels while 
top predators and basal species contribute little. Similar conclusions were further drawn from semi-analytical work 
j46j, where a maximum in biodiversity at an intermediate trophic level was predicted. Another consequence for 
such webs is that additions of species are generally not possible at any trophic level, if sustainable ecologies are to 
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be achieved. Even when the addition satisfies the assembly rules, its presence might cause substantial redistribution 
of biomass, i.e. shifts between consumer and resource limitation. In practice, it may be precisely these dramatic 
transitions that explain the profound and abrupt impacts on species abundance and energy flow patterns which 
are sometimes observed in the field. E.g. the introduction of opossum shrimp into a lake caused a cascade of 
trophic disruptions by reduction of salmon numbers and subsequent depletion of eagle and grizzly bear 147). On 
the other hand, our rules also describe the circumstances, under which removal of a species must trigger additional 
extinctions. 

The assembly rules thus allow predictions of secondary extinction, resulting either from addition or removal 
of species. If the modified food web obeys the assembly rules, the food web might be stable. Indeed, in some 
observed cases, ecological release of new species into a habitat has had relatively gentle effects 11481 . However, 
a violation of the assembly rules (Eas 1516b by addition of a new species can have one of two effects: Either the 
new species will not be competitive and collapse, or a number of species will collapse (possibly including the 
species itself) to restore the food web to a permitted state. For removal of a species that leads to violation of the 
assembly rules, secondary extinctions ll49l must be triggered to re-gain a sustainable state. We find a consistent 
pattern, when considering species removal in empirical food webs: E.g. in consumer-limited webs, such as the 
free-living empirical webs, secondary extinctions are more likely triggered by removal of resource than consumer 
species (Sec. S9 and Fig. S12). 

Community omnivory |j50] j2i :51! 521 and parasitism |[53l [27] |TT1 have been suggested as contributing to food 
web stability. Our approach provides theoretical support for this claim. We indeed find that full rank of the food 
web interaction matrix is difficult to achieve for species that are organized at strict trophic positions. Species 
that consume resources at different trophic positions, e.g. omnivores or some parasites, are shown to loosen the 
constraints and make it easier to comply with the assembly rules, i.e. finding a non-overlapping pairing of species. 
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Trophic Levels 



competitive exclusion possible coexistence not fulfilled 


Figure 1: Food web interaction matrix and application of perfect matching, a. White (gray) boxes indicate 
nonzero (zero) matrix elements, orange boxes are unity matrix elements for the primary producers; dark and light 
blue squares indicate a possible path chosen, allowing det(7l) to be nonzero. Here, N 0 = ri \ + 713 + W 5 = 18 and 
N e = ri 2 + 714 + Uq = 12, and Eq.[ 6 ]is fulfilled with A = 1. Inset: Schematic of a possible pairing for the chosen 
path. Note that the invariance property of det(7l) was used, yielding only m non-vanishing matrix elements in the 
lower right block (Details: SI), b. Perfect matching j28l applied to simple food webs where competitive exclusion 
rules out coexistence due to lack of niches (i) and where enough niches are available for coexistence (ii). (iii) and 
(iv) are two additional examples, where coexistence is ruled out by the assembly rules. In (iii), ni=2, 712 = 713 = 1 . In 
(iv), 77 , 1 = 713 = 714 = 1 , 712 = 2 . In both, A = N a — N e £ {0,1}, see Eq.0 
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Figure 2: Species richness for different trophic levels, a. Food web with two trophic levels only; a staircase 
of coexistence with balanced species richness at levels 1 and 2 |23l . b. Three trophic levels. The number of 
intermediate species must equal the total number of basal and predator species. Intermediate species dominate 
ecosystem biodiversity, c. Four trophic levels. 712 (W 3 ) must at least match basal (predator) species richness n\ 
( 714 ), indicated by thin black lines in green and orange bar. Solid green (orange) bars show the minimal upper 
bound to species richness in trophic level one (two). Species richness ri 2 and 713 can increase even further by 
co-evolution of intermediate species (shaded region). Note the applicable assembly rules shown for the different 
cases. 
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Figure 3: Food web assembly, a. Consumer limited food web and its interaction matrix. The symbols “x” mark 
nonzero entries and circles a path through the matrix. Basic nutrient is shown as a gray hexagon, whereas species 
on subsequent trophic levels are shown with coloured circles. Small circles highlight limitations by consumers. 
Shaded ovals indicate possible pairing of species, b. Possible assembly of the food web in fa) with labels for 
N a +N e and N a —N e . We set all growth and coupling constants equal to unity, but consider fine-tuned decay 
coefficients a<Cl J39). Note the transitions between biomass limiting states. Sizes of circles indicate approximate 
densities of species in the different states. Link shown in gray was set to zero in the numerical simulations in (c). c, 
Species added one-by-one as shown in (b) and granted small initial population density (10~ 4 ). After each addition 
the system is integrated until steady state is reached. In the plot, colors of curves denote species of similar colors 
in the panels of (b). Note the double-logarithmic axis-scaling in (c). Time in each panel is relative to the time of 
introduction of the new species. 
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Figure 4: Inclusion of omnivory. Food web and corresponding interaction matrix for several species with sharp 
trophic levels and a single generalist omnivore. We only show its paired link in the network plot, but indicate all 
its interactions in the matrix. 


24 


































: Fii!m;^;i5;ii;i ;'i 

->^-3s4v?Si! p-^-l r 

TfsaJS,. ''IJf :•■ 

[\ 


WL&L=o) 



■!*S3K3w®3£»g3» 

2 t | nr - rr ? <{ >-•-£.---: ** 

= vas&££jKS)i£i 

aamiUMJtlX 

irimg '•£&&*$■ 

s rfiisiiMMW^ l, 

= .mrr^rae^! 



ij-ljPS? 

M.V*! .... 

:iii 

If mill!!, 

A 1 



*■■'• i.'-;V k Yw> '*.1\:t;f* 

• .r y- • - ••*' ...•*..•• 



Figure 5: Comparison to existing models and data analysis, a, Bahia Falsa free-living food web, with blue and 
red matrix elements for predator respectively prey dependency. Trophic levels indicated by gray and black bars, 
whereas color coding along the left and upper edge labels chain length of free-living species (increase from red to 
blue shades), b. Niche model simulation of the Bahia Falsa free-living food web. c. Cascade model simulation of 
the Bahia Falsa free-living food web. d, As (a) but with parasites (“Par”) and with colors for free-living species 
carried over from (a). Remaining species are parasites, e. As (d) but additionally including symmetric concomitant 
links (“ParCon sym”). f. Lack of niches (d) for seven empirical food webs lUTI . Labels mark the sub-webs of free- 
living species (“Free"), including also parasite links (”Par“), asymmetric (’’ParCon asym“) and several symmetric 
concomitant links (’’ParCon sym“). Dashed line connects averages in these categories. Rank deficiencies for free 
living (df), free-living and parasite species (d/+ p ) as well as additional concomitant links (d CtSV m) marked in (a), 
(d), and (e), respectively. (Analysis details and abbreviations: see Methods). 
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Figure 6: Simulations of different food web matrices, a, Barplots indicate distributions of node richness for 
each approximate trophic level in the seven empirical foodwebs and a generic foodweb derived by averaging the 
empirical node richnesses in each trophic level, b. Interaction matrix corresponding to the generic food web, 
containing 110 free-living and 47 parasite species (Details: Methods and Sec. S10). c. Addition of parasites that 
form random links to any existing free-living species, d. Addition of parasites that are confined to consumer at a 
specific trophic level, e. Similar to (c) but with the restriction of parasites consuming only free-living species at 
levels 3 and 4. f. Similar to (e) but with additional parasite-parasite interactions (hyperparasitism), approximately 
5 percent of parasite links are from parasite to parasite. Note the color coding along the edges of the matrices in 
(c)—(f), chosen as in Fig.[5ji,d,e. g. The lack of niches, i.e. rank deficiency, as a function of the number of links 
per parasite for each of the four cases described in (c)—(f). 
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